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Abstract. We present an analysis of the excitation conditions of the molecular gas in the streamers and the 
outflow of M82 based on observations obtained at the IRAM 30 m telescope. Our analysis of J=1^0 and J=2^1 
transitions of CO and 13 CO and the CO(J=3— >2) line in 13 regions outside the central starburst disk shows that 
the gas density within the streamer/outflow system is about an order of magnitude lower (log(nn 2 ) ~ 3.0cm -3 ) 
than in the central molecular disk. We have used an LVG model and data from the literature to constrain the 
flux density in each CO transition (the 'CO line SED') arising from the streamer/outflow system and the central 
starburst disk itself. Globally, we find that the CO flux density up to the J=3— >2 line is dominated by the diffuse 
outer regions while lines above the J=5— >4 transition are almost exclusively emitted by the central starburst 
disk. We compare the CO line SED of M 82 to CO observations of galaxies at high redshift and suggest that 
small high-J/low-J CO flux density ratios (observed in some of these sources) are not necessarily caused by a 
different excitation of the central molecular gas concentration, but may result from an additional, more extended 
and diffuse gas reservoir around these systems, reminiscent of the situation in M 82. 
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1. Introduction 

Studying the physical properties of the molecular gas in 
galaxies is of fundamental importance for understanding 
\ the processes leading to star formation at low and high 
redshifts. Due to its proximity and its strong emission in 
molecular lines, the nearby starburst galaxy M 82 is one of 
the best studied starburst environments. The central con- 
centration of molecular gas, which feeds the strong star 
formation activity, has been the subject of many studies 
addressing both the excitation conditions and the distri- 
bution of the molecular gas (e.g., Young & Scoville H 9841 
Wild etal. 113321 Weifietal. [TUM Mao et al. I2UUU1 WeiB 
et al. 120011 Ward et al. 2003|. Some of these studies already 
have provided evidence for extended CO emission around 
M82's center (e.g. Stark & Carlson llTjgl Seaquist & Clark 
EUoTl Taylor, Walter, & Yun l2UUTjl . A high-resolution, 
wide-field CO study revealed that large amounts of molec- 
ular gas are present in molecular streamers and the out- 
flow around M82's central disk (Walter, Weifi & Scoville 

mm. 
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Recent studies of CO line and dust continuum emis- 
sion in distant quasar host galaxies (z > 2) have used 
M82's central starburst disk as a nearby template (e.g. 
Yun et al. 127177711 Cox et al. 127)011 Carilli et al. l2()()2al Carilli 
et al. I2002bl Weifi et al. I2003fl . At cosmological distances, 
however, only integrated values can be derived as usually 
no information on the spatial distribution of the molecu- 
lar gas exists. Therefore CO line ratios can only provide 
information on the galaxies average, global excitation con- 
ditions. Consequently a comparison of such line ratios to 
values determined only for the center of M 82 are incom- 
plete and biased as they neglect the contribution of a pos- 
sible extended molecular gas component which could po- 
tentially have significantly different excitation conditions. 
This is the main motivation for the present study in which 
we present an analysis of the excitation conditions of the 
extended molecular gas surrounding M82's central star- 
burst disk. 

2. Observations and data reduction 

Using the IRAM 30 m telescope we observed the 
CO(J=1->0), CO(J=2->l), 13 CO(J=1^0) and 
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Fig. 1. Logarithmic representation of the integrated CO(J=1^0) map of the zero-spacing-corrected OVRO mosaic 
(Walter, WeiB & Scoville 2002). Circles indicate the observed positions for which we analyzed the CO line brightness 
temperature ratios. The size of the circles corresponds to 22", i.e. the spatial resolution of the 30 m telescope at 3 mm. 
The white dashed contour shows our division between the central molecular disk and the streamer/outflow system. 



13 CO(J=2^1) transitions towards 13 positions cov- 
ering the molecular outflow and the molecular streamers 
surrounding M82's starburst toroid (see Fig. Spectra 
were obtained using the wobbler switch technique with a 
wobbler throw of 240" and a wobbler frequency of 0.5 Hz. 
The receiver alignment was determined using pointing 
scans on Saturn and was found to be better than 2". 
Pointing was checked every 2h on the nearby quasars 
0954+658 and 0836+710. We estimate the pointing 
accuracy to be better than 3". Calibration was checked 
once per observing run on W30H and the south-western 
molecular lobe in M82. 

System temperatures were on average T SJ/S «180K and 
330 K at 115 and 230 GHz respectively. The data were 
recorded using the Vespa auto correlator with 512 MHz 
bandwidth and 1.25 MHz resolution at 3 mm and two 
1MHz filter banks (512 channels each) at 1mm. The 
velocity coverage at 115 GHz is 1330 kms -1 with 3.4 
kms -1 resolution. The corresponding values at 230 GHz 
are 660 kms -1 and 1.3 kms -1 . The telescope's resolution 
at 230 and 115 GHz is 11" and 22" respectively. To enable 
comparison between the brightness temperatures at 
230 and 115 GHz we observed a cross consisting of five 
pointings at each position (separation 6"). Depending on 
the CO line strength the integration time at each central 



pointing was between 12 and 54min. All other pointings 
within the cross were observed for 6 or 12min. 
For data processing we used the CLASS software package. 
Only linear baseline were removed from each spectrum. 
For the CO(J=1->0) and the 13 CO(J=1^0) lines only 
spectra from the central pointing were coadded. For the 
reduction of the CO(J=2-^l) and the 13 CO(J=2^1) 
data, spectra at each pointing in the cross were coadded. 
From the resulting 5 spectra we generated the final 
spectrum at 22" resolution by a weighted addition with 
a weight of 1 for the central and 0.71 for the other four 
spectra (Johannsson etal. I1994[l . The coadded spectra 
were converted to the main beam brightness temperature 
scale (T mb = F cS /B cS T* A ) using F cff = 0.95, F cff = 0.91 
and B cS = 0.71, B cS = 0.49 at 115 and 230 GHz respec- 
tively (IRAM newsletter 55, 2003). Finally, CO(J=1^0) 
and CO(J=2-^l) spectra were smoothed to 3 kms -1 
resolution. The final resolution of the 13 CO(J=1^0) 
and 13 CO(J=2^1) spectra is 10 kms" 1 . In addition 
to theses observations we used the CO(J=3^2) data 
from Dumke etal. I|2001[) . To allow a comparison to our 
data we smoothed the CO(J=3^2) data cube to 22" 
resolution before extracting spectra at each position 
covered by our observations. Individual spectra at 22" 
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Fig. 2. Spectra of the CO(J=1->0), C0(J=2->1), and CO(J=3^2) (left), and the 13 CO(J=1^0) and 13 CO( J=2->1) 
(right) transitions at 22" resolution at individual positions. Line brightness temperatures are given in Kelvin on a T m b 
scale. The CO(J=3— >2) spectra have been calculated from the data in Dumke etal. 2001 



resolution are shown in Fig. [2J Table ^ summarizes the 
line parameters at each position. 

3. Results 

3.1. Line brightness temperature ratios 

The CO(J=1->0), CO(J=2^1) and 13 CO(J=1->0) lines 
are detected towards all positions. At position 13 
(streamer S4) the 13 CO(J=2^1) line is not detected (So- 
upper limit: Ti3 CO(J=2 ^ 1) < 20 mK). The CO(J=3^2) 



line is detected at all positions except positions 13 which 
is outside the field observed by Dumke etal. I|2001|l . The 
line profiles of all transitions at individual positions are 
in good agreement except for the northern outflow (Pos. 
4 & 5) where the wings of the CO( J=3— >2) line are broader 
than the corresponding profiles for the CO(J=l— >0) and 
CO(J=2^1) lines (see Fig. |2J. We attribute these line 
profile differences to small pointing offsets between the 
observations by Dumke etal. and our observations. 
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Fig. 2. continued. Note that the CO(J=3-^2) spectrum at position 11 has a offset of — 5",0" with respect to the other lines 
due the limited extend of the CO( J=3^2) cube. The spectra at the bottom right panels are average spectra from all positions. 



Table 1. CO line parameters at 22" resolution. 
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a CO(J=3^2) data from Dumke et al. I200T1 

b CO( J=3^2) line profile wider than CO( J=l— >0) and CO(J=2^1) line (see discussion in Sec. 13.11 
c CO( J=3^2) spectrum has a offset of -5", 0" from Pos. 11 
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We have used our high-resolution CO map (see Fig.^l 
and the beam patterns of the 30 m telescope (Greve, 
Kramer & Wild 19S5J) to estimate the contribution from 
the error-beam in each spectrum. Using the first two 
error-beam patterns we find that the contribution is neg- 
ligible at all positions (<5%). 

The integrated and peak brightness temperature ratios 
agree well within the observational errors. In the follow- 
ing we use the integrated brightness temperature ratios, 
since this provides a higher signal to noise ratio for the 
13 CO data. All R ratios used throughout this paper are 
defined in Tabled Only in the case of R31 in the north- 
ern outflow (Pos. 4 & 5) we used the peak brightness tem- 
perature ratio, as the integrated CO(J=3^2) intensities 
are too high due to the different line profile (see above). 
Statistical errors are negligible for all CO lines. From our 
W30H scans we estimate the systematic error of the Rio 
ratio to be « 15%. The different observing methods at 
3 mm and I mm lead to « 30% systematical errors for 
i?2i and i?2i- Taking the independent uncertainties of our 
CO(J=1^0) intensities and the CO(J=3^2) errors re- 
ported by Dumke et al. I|2001[l lead to « 30% error for the 
i?3i ratio. Line brightness temperature ratios and their 
errors, taking the statistical error for the 13 CO lines into 
account, are summarized in Table |2J 

Fig.|31shows the spatial variations of the most relevant 
line brightness temperature ratios. Differences between in- 
dividual positions are remarkably small. The only clear 
trend visible in the spatial distribution of the line ratios is 
a decreasing R31 with galactocentric distance. This trend 
is visible along the streamer SI and S3 as well as along the 
northern and southern outflow. Along S2, however, R31 is 
constant. All other ratios do not show a clear trend with 
increasing distance from the central starburst region. R21 
is around unity for all positions covering the streamer sys- 
tem, slightly above unity in the southern, but only about 
0.8 in the northern outflow. The corresponding ratios in 
13 CO, i?2 3 , show a larger scatter and values range between 
0.7 and 1.5. Judging from the ratio of the ground transi- 
tions of CO and 13 CO, Rio, the analyzed positions can be 
subdivided into positions within and outside the optical 
disk. Within the disk, along the molecular streamers SI & 
S2 and at the base of the northern and southern outflow 
the ratio Rio is in the range of 20 — 26. Outside the opti- 
cal disk, at the end the outflows and within the streamer 
system S3 and S4 the ratio rises to value between 29 — 39. 
While the line ratios at the base of the southern outflow 
almost resemble those in the molecular disk of M 82 (see 
Sec. 14.2(1 . the gas in the northern outflow appears to be 
much less excited. In fact, R21 and -R31 in the northern 
outflow are even lower than in the molecular streamers. 

3.2. Excitation Conditions 
3.2.1. Average Properties 

To analyze the excitation conditions of the molecular gas 
surrounding M82's starburst region we use a spherical, 



isothermal one component large velocity gradient (LVG) 
model. Given the small variations of the line brightness 
temperature ratios with position we use in the following 
the line ratios determined from the average spectra to in- 
vestigate the global excitation away from the star forming 
regions. Note that these ratios differ from the numerical 
average of the line ratios at individual positions since they 
correspond to the intensity weighted averages which would 
be observed if M 82 was shifted to cosmological distances 
and the emission would remain unresolved. Spatial varia- 
tions of the excitation conditions are discussed in section 
13.2.21 The averaged spectra are shown at the bottom right 
panels of Fig. [3 the corresponding averaged integrated 
line ratios are given in the bottom row of Table El 

We compare LVG predicted line ratios to the obser- 
vations in a LVG parameter space covering densities and 
kinetic temperatures of log(nn 2 ) = 2.0 — 5.0 cm~ 3 and 
?kin = 5 — 200 K. CO abundances per velocity gradient 
and 13 CO abundances relative to CO are in the range 
of [CO]/dv/dr = 1 • 10~ 6 - 2 • 10~ 4 pc ( kms" 1 )- 1 and 
[CO]/[ 13 CO] = 30-100. 

The observed CO ratios, -R21 and R31, limit solutions 
in the LVG parameter space to a density range between 
log(nn 2 ) — 2.3 — 3.5cm -3 where the lower density limit 
corresponds to the highest value of [CO]/dv/dr and the 
upper density limit to the lowest [CO]/dv/dr. The mea- 
sured high Rio ratio restricts solutions to physical condi- 
tions with moderate opacities of the CO(J=1^0) (T12 ~ 
1 — 3) and very low opacities in the 13 CO( J—l— >0) tran- 
sition (n3 « 0.01 — 0.05) within this density range. As 
a consequence solutions are limited to small CO abun- 
dances per velocity gradient (< 2 • 10 -5 pc ( kms -1 ) -1 ) 
and [CO]/[ 13 CO] abundance ratios > 50. This additional 
constraint on [CO]/dv/dr yields an allowed density range 
of log(nn 2 ) = 2.7 — 3.5 cm~ 3 and small CO column den- 
sities per velocity interval (< 2 • 10 17 cm~ 2 (kins -1 ) -1 ). 
The kinetic temperature is not well constrained without 
assumptions on the values of the [CO]/[ 13 CO] abundance 
ratio and the CO abundance per velocity gradient. A lower 
limit on the kinetic temperature of 7ki n > 30 K is given by 
the densest possible solution within the parameter space. 

In Figs .01 & we show LVG solutions for the aver- 
aged line ratios as a function of abundance and velocity 
gradient variations. The lines represent the temperature- 
density combinations for which the LVG predicted line 
ratios match the average observed line ratios. The fig- 
ures show that lower CO abundances per velocity gradient 
and/or higher fractional 13 CO abundances lead to cooler, 
denser solutions while higher abundances per velocity gra- 
dients and/or lower fractional 13 CO abundances lead to 
warmer and more diffuse solutions, respectively. For a CO 
abundance per velocity gradient of 1 • 10~ 5 pc ( kms -1 ) -1 
and a [CO]/[ 13 CO] abundance ratio of 80 the averaged line 
intensity ratios are well reproduced by a kinetic temperate 
of Tkin = 55 K with a H2 density of log(nn 2 ) = 3.2 cm -3 . 
These values correspond to the intersection of the lines 
shown in Fig.0]&|3] (center). 
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Table 2. Selected integrated line brightness temperature ratios at 22" resolution. Errors include 15% systematic 
uncertainty for R w , 30% for i? 2 i and R^f and 30% for i? 32 . Statistical errors have been included for the 13 CO lines 
only. 
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a Peak brightness temperature ratio (see text). 

b CO( J— 3— »2) spectrum has a offset of —5", 0" from Pos. 11 

c Intensity weighted average determined from the integrated line brightness temperature ratios of the averaged spectra (see 
Fig-IH bottom right panels). 
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Fig. 3. Spatial variations of the R21 (upper left number), R31 (lower left), R 21 (upper right) and i?io (lower right) 
line brightness temperature ratios (see Table [21 for definitions). The greyscale is an optical V-band image of M82 
(Cheng etal. 1997). The contours represent the distribution of the CO(J=1^0) emission (see Fig. P). 
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Fig. 4. LVG line brightness temperature ratios as a function of the kinetic temperature and H2 density for a fractional 
13 CO abundance of 60, 80 and 100. The CO abundance per velocity gradient is 1 • 10 -5 pc (kms -1 ) -1 for each graph. 
Displayed line ratios correspond to the averaged observed values listed in the bottom row of Table 
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Fig. 5. LVG line brightness temperature ratios as a function of the kinetic temperature and H2 density for a CO 
abundance per velocity gradient of 2 • 10 -5 , 1 • 10" 5 and 5 • 10~ 6 pc (kms -1 ) -1 . The fractional 13 CO abundance is 80 
for each graph. Displayed line ratios correspond to the averaged observed values listed in the bottom row of Table El 



3.2.2. Individual regions 

In the following we compare the excitation conditions 
at individual positions. To do so, we assume that the 
abundance and velocity gradient variations do not change 
with position and fix the CO abundances per veloc- 
ity gradient and the fractional 13 CO abundance to 1 • 
10 -5 pc ( kms" 1 ) -1 and 80 respectively. 

Along the streamers SI & S2 the excitation conditions 
are fairly similar. The decrease of the R31 along SI is re- 
flected in a small decrease of the density and temperature 
from log(nn 2 ) = 3.2cm -3 and T^ n — 50 K at the base 
of the streamer to log(nn 2 ) — 3.1cm -3 and Tki n = 40 K 
at its end. Similar conditions are met along streamer 
S2 {log(n H2 ) « 3.3 - 3.1cm -3 and T kin « 45 - 40 K). 



The base of the northern outflow (Pos.4) can be fit well 
with cooler gas (Tk; n ~ 30 K) at comparable density 
(log(nn 2 ) — 3.1cm -3 ) while the high Rio ratio com- 
bined with the low CO ratios, R21 and R31, at the end 
of the outflow leads to T^ n w 40K with a low density 
l°g{ n H 2 ) — 2.8cm -3 . In contrast, the base of the south- 
ern outflow consists of denser gas at higher temperature 
(Tkin ~ 70 K, log(nn 2 ) rs 3.5 cm -3 ), comparable to condi- 
tions in the starburst disk (see Sec. l4.2f> . At the end of the 
southern outflow as well as along the streamers S3 and S4 
(the positions with the highest Riq ratio) a warm diffuse 
medium is required to match the observed line ratios. The 
temperatures and densities range from Tki n w 60 — 120 K 
and log(n\i 2 ) = 2.7 — 3.0cm -3 respectively. Note that the 
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* 2 CO(2-1)/' 2 CO(1-0) ,2 CO(3-2)/ 12 CO(1-0) 

,2 C0(1 -0)/ 13 CO(1 -0) ,2 C0(2-1 )/ 13 C0(1 -0) 

— - ' 2 CO(2-1 )/ ,3 CO(2-1 ) 




2 2.5 3 3.5 4 



log[n(H 2 )] [cm 3 ] 

Fig. 6. LVG predicted R 2 i, R31, Rio, Rl 2 / 13 and i?^ 13 
line ratios as a function of H2 density and kinetic tem- 
perature for a CO abundance per velocity gradient of 
1 • 10 -5 pc ( kms -1 ) -1 and a fractional 13 CO abundance of 
80. Numbers indicate the respective line ratios. The sym- 
bols show the density - temperature combination of the 
streamers SI & S2 (grey circle), the northern outflow (grey 
square) , the base of the southern outflow (white square) 
and regions south of the optical disk (S3, S4, and the end 
of the southern outflow; grey ellipse). 

large difference in R21 and R31 at the end of the southern 
outflow (Pos. 7; 1.1 and 0.5 respectively) can not be fit by 
a single gas component. In Fig.Elwe show where the dif- 
ferent regions lie in the LVG temperature density space. 
Not all line ratios at a given positions can be fit equally 
well and the locations of the regions show roughly the best 
agreement with the measured line ratios. 

4. Discussion 

4.1. Properties of the streamer/outflow gas 

The high isotope line ratio observed in the J = 1 — 
transition (R\o) at all positions clearly rules out that 
the gas conditions in the outer parts of M82 are simi- 
lar to those in the disk of the Milky Way or other spiral 
galaxies. Typically the CO( J=l— >0) line has large optical 
depth in disk clouds in spiral galaxies and even the emis- 
sion in 13 CO(J=l— >0) is only moderately optically thin 
which results in Rio typically about 5. R\o much closer 
to the [CO]/[ 13 CO] abundance ratio is a well studied sig- 
nature of molecular gas in the centers of active galaxies 
(including M 82's starburst disk with Rio ~ 20) and merg- 
ers (e.g. Aalto et al. ITMT1 Sage & Isbell ITMT1 Henkel & 



Mauersberger 1199311 . The reason for the 13 CO depression 
in luminous starburst galaxies is still under debate but 
most explanations involve 13 CO/CO abundance ratio vari- 
ations and opacity effects (see e.g. Henkel & Mauersberger 
ITMfl Aalto et al. IT5§5l Kikumoto et al. IT53SI Tanigucchi, 
Ohyama & Sanders [TUM Paglione et al. I2UTT1 Meier & 
Turner 2004). Given its origin from the star forming disk 
and the likely interaction with outflowing ionized gas, high 
Rio ratios are therefore expected from the gas located in 
northern and southern molecular outflow. The observed 
increase of Rio along both outflows supports the idea that 
'superwinds' have large isotope line ratios ( Tanigucchi, 
Ohyama & Sanders 1999). For gas along the streamers SI 
and S2 the high R w ratios are unexpected as Rio is usually 
found to decrease with galactocentric distance (Paglione 
et al. I2001|l and since both streamers are located well out- 
side the central star forming regions. We note, however, 
that decreases of Rio typically only becomes significant for 
galactocentric distances larger than 2 kpc (Paglione et al. 
2001). Increasing Rio ratios in the vicinity of starburst 
have also been observed in e.g. IC 342 (Wright et al.[l993) 
and NGC6946 (Meier k Turner l2lTOjl . In these cases the 
large isotopic line ratio has been attribute to dispersed 
low-density gas in the 'spray regions' of barred potentials. 
Given the low density in SI & S2 and that these streamers 
form the base of the disrupted Hi features within M82's 
disk (Yun et al. I1993|l we speculate that their large Rio 
may result from dispersed gas due to tidal forces along SI 
& S2. 

In the context of our LVG models, the large isotopic 
line ratio is mainly a result of the very low optical depth of 
the ground transition of 13 CO (n 3 » 0.01 - 0.05). These 
low opacities result from the low 13 CO abundances per 
velocity gradient of our models. Outside the optical disk, 
in regions with Rio > 30, higher kinetic temperatures and 
lower H2 densities (resulting in lower column densities per 
velocity interval) push 13 CO opacities to the lower end 
of the range. For gas in local thermodynamic equilibrium 
(LTE) and large optical depth in CO(J=1^0) the ratio is 
given by Rio ~ t^ 1 . Although (i) the excitation temper- 
atures between CO and 13 CO differ up to a factor of 2 in 
our models (non LTE) and (ii) the moderate opacities in 
CO(J=1^0) introduce an additional deviation from the 
above equation, the opacity in 13 CO is the most impor- 
tant quantity for variations in Rio. This suggests that the 
CO emission outside the optical disk (streamers S3 & S4 
and at the end of both outflows) arises from diffuse and 
presumably warm gas with low opacities in 13 CO. 

Interestingly, the southern and northern outflows have 
different excitation conditions (mainly reflected in the low 
R21 and i?3i ratios in the northern outflow). The lower 
temperatures towards the north might indicate that the 
gas is less affected by outflowing ionized gas. Indeed, the 
outflow seen in Ha and X-rays is not symmetrical with 
respect to the optical disk but brighter and also more 
extended towards the south (e.g. Lehnert, Heckman & 
Weaver 1999). This, however, might simply be an effect 
of the orientation of the outflow with respect to the ob- 
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server since the southern outflow is inclined towards us 
(see e.g. McKeith et al. Greve 111)551 for the outflow 

geometry) . 

Although the kinetic temperatures of the molecular 
gas are difficult to constrain without better estimates of 
the CO abundances, the general picture to emerge from 
the LVG models is that most of the gas in the outer re- 
gions of M 82 has low densities and low optical depth at 
temperatures in excess of 30 K. 



4.2. Comparison to M82's center 

The physical conditions of the gas within the central 
molecular toroid have been the subject of many stud- 
ies (e.g. Wild et al. Glisten et al. HMfl Mao et al. 
I2T)0T)I Petitpas & Wilson I27TO1 Weifi et al.ESHI Ward et al. 
2003). These studies have shown, that the observed line 
ratios can be reproduced by emission from a low (LE) and 
a high excitation (HE) gas component (Glisten et al ll993l 
Mao etal. HUM Weifi etal. IMiTl Ward etal. HUU^i . In 
this picture the high excitation component, responsible for 
the excitation of levels beyond the J = 4 — > 3 transition, 
arise from dense gas at temperatures similar to or higher 
than the dust temperature (log(nn 2 ) w 3.5 — 4.5 cm~ 3 ; 
2kin > 40 K) . The low excitation component is emitted by 
diffuse gas at much lower density {log(nn 2 ) ss 3.0 cm -3 ). 

Our analysis shows that the gas surrounding M82's 
starburst region has similar properties as the low exci- 
tation component in the starburst center itself. Similar 
conclusions have been reached from an analysis of the CO 
J=3->2/J=2->l, 13 CO/CO J=3^2 and C 18 Q/ CO ra - 
tios in the molecular "halo" by Seaquist & Clark (2001). 
They interpret their decreasing line ratios as decreasing 
CO excitation and optical depth with increasing distance 
from the nuclear region. These findings are in line with 
our results. Seaquist & Clark find that the variations of 
the observed line ratios in the transition region between 
the center and the "halo" are consistent with a decreas- 
ing filling factor of the high excitation component with 
distance from the center. In this picture the CO emission 
from the outer regions arise exclusively from gas with sim- 
ilar properties than the low excitation component in the 
center. We discuss this scenario in more detail in the next 
section. This suggests that mainly the diffuse low density 
molecular gas in the starburst region is involved in the 
outflow. Given the higher inertia of dense (and therefore 
more massive) gas clouds, we speculate that the superwind 
affects the kinematics of gas at low density more strongly. 
However, the high CO excitation found at the base of the 
southern outflow indicates that also gas at higher density 
might be pushed out of the starburst regions, at least to 
scale heights of a few hundred pc, less than half the dis- 
tance seen in the outflow of the diffuse molecular gas. 



4.3. CO-Line-SED 

4.3.1. High J CO fluxes 

We have used published maps in various CO transi- 
tions towards M82 to estimate the total flux density 
emitted in each line from the center and the outer 
parts of the CO distribution. Large-scale maps (cov- 
ering the streamer/outflow system) only exist up to 
the CO(J=3— >2) transition. The observed flux densities 
from the outer parts of M 82 were determined by mask- 
ing the central molecular disk in the CO(J=1^0) and 
CO(J=2^1) IRAM 30m data cubes from Weifi etal. 
p)01|l and the CO( J=3^2) data cube from Dumke et al. 
H2()()l|l (after smoothing all data cubes to 22" spatial reso- 
lution and regridding onto the same grid). The mask was 
chosen to contain a flux density of w 5000 Jy kms -1 in 
the CO(J=1^0) data cube (see Fig^]), close to the value 
determined by us for the central molecular disk (Walter, 
Weifi & Scoville ED02) ■ CO flux densities from the center 
were determined within the inverted mask. For the center 
we also analyzed the CO( J=4^3) and CO( J=7— >6) cubes 
published by Mao etal. lj2lTOjl . Since the CO(J=4^3) 
map is larger than our mask we applied the same tech- 
nique as for the lower J lines to determine the central flux 
density in CO(J=4^3). The CO(J=7->6) map is smaller 
than the central mask and we adopt here the total flux 
density in the data cube to be representative for the flux 
density from the center. The central flux density of the 
CO( J— 6-^5) transition was taken from the map by Ward 
etal. (2003), which roughly covers the central mask. The 
observed CO flux densities for the center and the outer 
parts are summarized in Table|3| 

4.3.2. CO line SED of the starburst disk 

In Fig.d (left) we show the observed CO line SED (CO 
integrated flux density vs. the rotational quantum num- 
ber) in the center of M82. CO flux densities rise up 
to the J=6^5 transition, although we can not rule out 
that the true CO(J=7^6) flux density is higher due to 
the limited coverage of the J=7— >6 map. In this plot 
we also show the LVG predicted CO line SED for a 
high (HE) and a low excitation (LE) component in the 
center. Following the approach by Seaquist & Clark we 
used the LVG model derived from the average line ra- 
tios in the steamer/outflow system as LE component (LE: 
log(nn 2 ) = 3.2 cm~ 3 , T^m = 55 K, CO abundance per 
velocity gradient 1 • 10~ 5 pc ( kms -1 ) -1 ; see Sec. l3.2.f|) . 
Adding a HE component with log(nn 2 ) — 4.2cm -3 , 
Tkin = 50 K, similar to conditions in the molecular lobes 
(e.g. Weifi et al. 120010 . results in a good fit to the observed 
CO line SED of the center. For these models we adopted 
a filling factor of 1/20 and 1/15 for the HE and the LE 
component in a solid angel of r2 center « 1000 arcsec 2 (as 
defined by our mask for the center). For comparison we 
show the CO line SED in the center of M 82 predicted by 
a high resolution LVG analysis based on the J = 1 — and 
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Table 3. Observed and model predicted CO integrated flux densities in the center and outer parts of M82. 
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Ref. 1: Weifi et al. ISOuTl f30 m data), 2: Dumke et al. I2U0T1 3: Mao et al. I3U001 4: Ward et al. I30U31 

a total flux covered by the observations 



J = 2 - 1 transitions only (Weifi et al. I20U1|I . The fluxes 
have been calculated taking all individual LVG solutions 
within the center into account. This example shows, that 
LVG models based on the low-J lines of CO and their 
rare isotopes give surprisingly good estimates of the flux 
density in the high-J CO lines. 

4.3.3. CO line SED of total emission 

In a similar way we have calculated the CO line SED for 
the streamer/outflow system. Here we used the LE com- 
ponent with an area filling factor of 1 /30 (derived from the 
average CO( J=1^0) spectrum at 22" resolution and the 
LVG predicted line brightness temperature) and a source 
solid angle of Pouter ~ 14000 arcsec 2 (estimated from the 
zero-spacing corrected CO(J=l— >0) OVRO map). The 
predicted line flux densities are shown together with the 
observed values in Fig.[7] (right). The plot demonstrates 
that this simple two phase model for the center and the 
outer regions is in good agreement with the observations. 
More detailed studies of the gas in the center show that the 
LE component is most likely warmer than the 55 K used 
here (e.g. Weifi et al. lSUUTl Ward et al. l3UU3|) . Higher tem- 
peratures at sightly lower density, however, do not change 
the predicted flux densities in the LE component signifi- 
cantly. In any case the line SED in the center is dominated 
by the dense gas that excites the high- J CO lines. This is 
in line with CO models for the center of NGC 253, where 
the CO line SED is fit by a single high excitation gas com- 
ponent (Bradford et al. ESUUl Bayet et al. l2"UU3|l . 

The combined CO line SED from the center and the 
streamer/outflow system, which would be observed if 
M82 was shifted to cosmological distances, however dif- 
fers strongly from that derived from the center alone (see 
Fig. right). The large area of the surrounding gas makes 
it the dominant component in the emission of the low- 
J CO transitions. Since this gas is not dense enough to 
produce significant flux beyond the CO(J=5^4) line, its 
CO line SED peaks already at the CO( J=3^2) transition. 
Consequently, even though two-thirds of the CO( J=l— >0) 



Table 4. LVG model predicted CO fluxes in the center 
(fi = 2.35 • 10~ 8 sr) and outer (O = 3.29 • 10~ 7 sr) parts of 
M82. 
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CO(2-l) 


1.9 


5.1 


6.9 


CO(3-2) 


5.3 
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17.2 


CO (4-3) 


9.8 


13.9 


23.8 


CO(5-4) 


15.0 


8.7 


23.7 


CO(6-5) 


20.2 


3.1 


23.3 


CO(7-6) 


21.2 


0.8 


22.0 


CO (8-7) 


14.0 


0.1 


14.2 


CO(9-8) 


4.6 


0.01 


4.6 



flux density arises from regions outside the star forming 
disk, CO lines above CO(J=5— >4) transition are mainly 
emitted from the center. In the combined CO line SED, 
flux densities rise up to the J=4^3 transition and de- 
crease beyond this line. We summarize the CO flux densi- 
ties derived from the LVG models in Table [21 Table |3] also 
gives the fraction of the total flux density emitted by the 
center in each CO transition. For comparison with other 
studies we give the model predicted line fluxes in Table 0] 

4.4. CO line SEDs at high redshift 

Currently CO line SEDs in sources at redshift > 2 are 
still poorly constrained. A few sources, however, have 
been observed in more than two CO transitions which 
allow at least for a quantitative comparison with the 
CO line SED of M82. Among the best studied exam- 
ples are the Cloverleaf quasar (z = 2.5, Barvainis etal. 
ITW1 Weifi etal. l2T)T)31l . IRAS FSC 10214 (z = 2.3, 
Brown & Vanden Bout 1992, Solomon, Downes & Radford 
HSU, PSS 2322+1944 (z = 4.1, Cox et al. IMil Carilli 
etal l2002a|l, BR12 02-0725 (z = 4.7 Omont et al. H"9"9l 
Carilli et al. l2002b|l and J1148+5251 (z = 6.4, Walter et al. 
I2T)TI31 Bertoldi etal. [2003). In all these sources the CO 
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Rotational Quantum Number J Rotational Quantum Number J 

Fig. 7. Left: Observed and model predicted CO line SED for the center of M 82 using a two component LVG model. The 
dotted line shows the predicted CO line SEDs of a high (HE), the lower dashed line of a low (LE) excitation component. 
The LE model correspond to the physical conditions determined in the outer regions of the CO distribution, the HE 
component is similar to the conditions in the molecular lobes. The black line is the sum of both components. For 
comparison we also show the CO line SED for a multi component LVG model (dashed line, Weifi etal. [2001). The 
filled squares are those CO observations which have been used to constrain the multi component model. Open squares 
are CO flux densities compiled from the literature. Right: CO line SED for M82 including the molecular streamers 
and outflow. The dashed lines show the model for the center (same as solid line to the left) and for the streamers and 
outflow. The solid line shows the total predicted CO line SED in the central 3x3 kpc. Filled rectangles are observed 
flux densities from the center of M82, open circles the observed flux densities in the streamer /outflow system and 
filled circles the total observed flux density in the central 3x3 kpc of M82. Observed and model predicted line flux 
densities are given in Tabled! 



flux densities are rising with increasing rotational quan- 
tum numbers (at least up to the J = 6 — > 5 transitions), 
similar to the center of M82 and NGC253. These studies 
argue against large extended low excitation CO halos sur- 
rounding these objects at high redshifts. However, there 
are a few examples where the flux density from the CO 
J=3— >2 transition is similar or higher than that of the 
J=7^6 line (see e.g. SMM 14011+0252 (z = 2.5, Downes 
fc So lomon 12717731 . SMM 16368+4057 (z = 2.4, Neri etal. 
2003 ) . In light of our new results in M 82 we suggest that 
the low J=7^6/J=3^2 ratios in these sources may not 
necessarily be caused by a different excitation of the cen- 
tral molecular gas concentration, but may result from an 
additional more extended and diffuse gas reservoir. 

5. Summary 

We have presented observations of the J=1^0 and 
J=2^1 transition of CO and 13 CO towards 13 selected 
regions outside the central starburst disk of M 82 cover- 
ing the prominent molecular streamers and the CO out- 
flow. Our observations show that the J=2^1/J=1^0 
and J=3^2/J=1^0 ratios of CO are lower in all 



streamer/outflow regions than in the central molecular 
disk. The CO ratios are fairly constant along the ana- 
lyzed positions and show only a slight trend of decreasing 
ratios with increasing distance from the center. The iso- 
tope line ratio CO( J=1^0)/ 13 CO(J=1^0) is around 25 
along the most prominent molecular streamers, similar to 
values found in the starburst disk, and increases along 
the outflows to 35. Similarly high values for this ratio are 
found in the diffuse streamers S3 and S4 south-east of the 
optical disk of M82. 

An LVG analysis of the gas suggests that the low 
CO ratios combined with the high isotope line ratio are 
best explained by a low density (log(n}i 2 ) « 3.0 cm~ 3 ) 
molecular environment with kinetic temperatures in ex- 
cess of 30 K and [ 13 CO]/[CO] abundance ratios - 80. 
In this picture the high CO( J=1^0)/ 13 CO( J=l->0) ra- 
tio is mainly caused by very low opacities in the 13 CO 
lines. Our analysis does not allow to rule out denser 
gas at lower kinetic temperature if abundance variations 
are taken into account. However densities in excess of 
log(nn 2 ) ~ 3.5 cm~ 3 require [CO]/[ 13 CO] abundance ra- 
tios > 100 or very small CO abundances per velocity gra- 
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client ([CO]/dv/dr< 1 • 10~ 6 ) to explain the observed line 
ratios which we consider unlikely. 

We use the LVG predicted line brightness temperature 
ratios of the CO J >3 transitions to estimate the total flux 
density emitted in the high-J lines from the diffuse CO 
component surrounding M82 starburst region. We find 
that the density of the streamer/outflow gas is not high 
enough to produce significant flux above the J=5^4 line. 
As a consequence the CO-line SED of the outer gas drops 
off beyond J=3^2 line while it is rising at least up to 
the J=6^5 line in the central molecular disk. The total 
J=1^0 and J=2^1 flux density within the central 3kpc, 
however, is clearly dominated by the emission from the 
outer regions. This implies that one has to exercise cau- 
tion when comparing global, unresolved, high-z line ratios 
to the values measured in the centers only of nearby star- 
burst galaxies. We conclude that small high-J/low-J CO 
ratios ratios seen in some high-z sources are not necessar- 
ily caused by a different excitation of the central molecular 
gas concentration, but may result from an additional more 
extended and diffuse gas reservoir around these systems. 
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